On the large order asymptotics of the wave function 

perturbation theory 



O.Yu. Shvedov 

Institute for Nuclear Research of the Russian Academy of Sciences, 
60-th October Anniversary Prospect 7a, Moscow 117312, Russia 

and 

Moscow State University 
Vorobievy gory, Moscow 119899, Russia 

February 1, 2008 
Abstract 

The problem of finding the large order asymptotics for the eigenfunction perturba- 
tion theory in quantum mechanics is studied. The relation between the wave function 
argument x and the number of perturbation theory order k that allows us to construct 
the asymptotics by saddle-point technique is found: x/k l l 2 = const, k — > oo. Classi- 
cal eualidean solutions starting from the classical vacuum play an important role in 
constructing such asymptotics. The correspondence between the trajectory end and 
the parameter x/k 1 ^ 2 is found. The obtained results can be applied to the calculation 
of the main values of the observables depending on k in the k-th order of perturbation 
theory at larges k and, probably, to the multiparticle production problem. 





It is well-known that the perturbation series in quantum mechanics, as well as in 
quantum field theory, diverges JlJ because of the factorial growth of the coefficients at 
large orders. High order asymptotics of the coefficients, which is important, for exam- 
ple, for the problem of the divergent series summation, has been discussed in many 
papers. For the quantities such as the ground state energy in quantum mechanics 

|3| or Green functions in quantum field theory [|j], the asymptotics is determined 
by the euclidean finite action classical solution satisfying the following boundary con- 
ditions: when "euclidean time" r tends to ±00, the euclidean solution tends to the 
classical vacuum. 

In this paper the role of classical euclidean solutions with other boundary condi- 
tions is discussed. It occurs that when the large order asymptotics not for eigenvalue, 
but for eigenfunction in quantum mechanics is considered, one must treat the classical 
trajectory starting from the classical vacuum as r — > —00 but finishing at some point 
Q as r = 0. 

Namely, consider, as usual, the following dependence of the Hamiltonian on the 
perturbation theory parameter g, 



where A is the Laplace operator, x = (xi, ...,x n ) is a set of coordinates, function V 
is a potential, which has a local minimum as its argument is equal to zero. Without 
loss of generality, suppose that V(Q) ~ Q 2 /2 + 0(Q 3 ) as Q — > O.The dependence (D) 
was considered in [fj, || . 

Consider the perturbation theory for the ground state eigenfunction. This per- 
turbation theory can be treated in several ways, namely, one can change the variable 
gx = Q and consider the wave function asymptotic expansion at fixed Q; this expan- 
sion is a tunnel semiclassical expansion because of the dependence (|l|), the analog of 
the Planck constant is g 2 . Large order asymptotics of semiclassical expansion can be 
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also found ||, and the asymptotics of the ground state energy perturbation theory, 
which was obtained in |2|, |3|, ||, [7]] for various cases can be checked || . 

Nevertheless, consider in this paper the ordinary way of treating the eigenfunction 
perturbation theory, 

oo 

vEr(x) = Xy* fc (z). (2) 

fc=0 

The purpose is to find the ^ asymptotics at large k. In order to do it, let us write 
through the functional intrgral. Make use of the path integral representation of 
the ground state wave function, 

= J Dqexp ^-l%g]j , 

where the integral is taken over the trajectories q(r) satisfying the boundary condi- 
tions, q(— oo) = 0, g(0) = x, and S is the euclidean action of the theory, 

S[Q] = J dr[Q 2 /2 + V(Q)]. (3) 

One can develop the perturbation theory for the ground state wave function with the 
help of the expansion of the integrand in the path integral representation into a series 
in g, the k-th order of this expansion can be written as a contour integral, so the k-th 
order of the wave function is given by 

the contour C(q) runs around the origin counterclockwise. 

Let us transform the integral @ to a saddle-point form. For this purpose rescale 
the variables q(r) and g in \fk times: 

g = Vl/Vk, q{r) = £(r)Vk. (5) 

This substitution means that the boundary condition for q,q(0) = x, must be rescaled 
in \fk times, too. Thus, the ^ asymptotics at large k and fixed x is not of a saddle- 
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point type, while such asymptotics as 

k — > oo, x/Vk — > £ = const (6) 

is obtained from the rescaled integral with the help of saddle-point technique and is 
given by the following expression, 

* fe (£oVfc) ~ T(k/2) expi-kAfo)), (7) 

where the function A equaling to 

A(Zo) = ls[V\$ + ~(\n±-l) (8) 

is taken in the extremum point, £ is a trajectory starting from zero as r — ► — oo and 
finishing at the point (o as r = 0. Varying the trajectory £(t), one can obtain that 
the function Q(r) = £(t)vA satisfies the classical euclidean equation of motion, 

Q = W(Q). (9) 

When one varies the quantity A, the following equation for it can be obtained: 

A/2 = J dr[V(Q) - QVV(Q)/2], (10) 

where r G (— oo,0]. We can treat the parameter ( |10D as a function of the end 
2(0) = Q of the classical trajectory <2(r) satisfying eq.(§). So we obtain the following 
relations connecting the point Q and the parameter £ in eq.(^), 

Q/^KQ) = ^. (ii) 

When one searches for the asymptotic formula for ^ under conditions (^), one must 
first find the boundary condition for the classical euclidean solution at r = from 
eq. ([TT|) , then find the classical trajectory arriving at this euclidean time moment in 
point Q and, finally, substitute the classical euclidean action on the trajectory and 
the parameter A to eq.(^). Notice that the function A in eq. (|TT|) can be multivalued. 
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The correspondence (|TT|) is non-trivial and interesting. Namely, suppose the tra- 
jectory Q(t) to spend all euclidean time in the vicinity of the point Q = 0, so that the 
coordinates of the trajectory end are small. Then the parameter A is of order 0(Q 3 ). 
Therefore, the corresponding value of £o is large. Thus, if one makes an attempt 
to find the behaviour of the asymptotics at large £o> trajectories closed to the path 
Q = must be taken into account, so one is in need only of knowledge of the first 
non- vanishing term of the difference V — Q 2 /2 at small Q. This statement can be 
illustrate as follows. In the each order of the perturbation theory the wave function is 
expressed as a product of the 3/c-th order polynomial in x function by the Gaussian 
exponent. It is the term of order x 3k expressed only through the third derivatives 
8 3 V/dQidQjdQk of the potential V at Q = that becomes the main one at large £ - 

We have seen that classical euclidean solutions starting at r — > — oo from zero 
correspond to large order behaviour of while the parameter £o in the condition 
(H) depends on the coordinates of the trajectory end. We can consider all such points 
in Q-space that the discussed classical trajectory passes through them. According 
to eq.flTT|), each such point Q can be associated with the corresponding point in £ - 
space by shifting of the parameter r to a constant in order to make the euclidean 
time moment at which the trajectory passes through the point Q to be equal to zero. 
Therefore, each trajectory in Q-space is related with the trajectory in £ -space (fig.l), 
£o( r ) = Q(r)/y A(Q(r)). Focus the attention on the distinction between this path 
and the trajectory £(r) = Q(t) / \J X(Q(0)) . According to the previous discussion, the 
trajectory £o( r ) starts from infinity, while the trajectory £(r) starts from the origin. 
Another interesting feature of the path £q(t) is that the derivative of the function A 
with respect to r is written in the form, A(£q(t)) = 7r £o, where tto = Q/VX- The 
form is analogous to the formula dS = PdQ. 

If the trajectory in Q-space finishes as r — ► +oo at zero then the corresponding 
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curve in £ -space finishes in the point £o = as r — > +00, too, because the integral fllCf ) 
becomes the classical action at larges r, which is not equal to zero. Thus, when one 
makes an attempt to find the asymptotics at small £0; one must consider the classical 
trajectory closed to the loop-like euclidean solution (fig.l) with the action S . Since 
the function A is approximately equal to 2S up to 0(Q 3 ), the function A is equal to 
I In iSq — to the same accuracy. The ^k asymptotics as fc|£| 3 <C 1,A;|£| 2 ^> 1 is, 
therefore, given by T(k/2)S k l 2 e k( > 2 1 2 . Notice that the pre-exponential factor in the 
asymptotics, which is omitted in this formula, diverges. Namely, it is connected 
with the determinant of fluctuations near the classical solution. When Q is small 
and the trajectory is closed to the loop-like solution, we come across the "almost zero 
mode" of the solution: when we shift the parameter r to a large constant, Q(0) is 
shifted to the small quantity. So the determinant is small. Therefore, we can expect 
the following form of the ^k asymptotics at fixed x and large k, 

^xM. (12) 

where x( x ) is n °t square integrable because of its exponential growth at larges x, 
x(x) ~ e^ 2 / 2 . 

Since the ground state eigenvalue is defined by the behaviour of the eigenfunction 
near the point £0 = 0, large order asymptotics of perturbation theory for it has the 
form similar to ( p"2|) up to a pre-exponential factor. This formula coincides with the 
results obtained in [Q, |, |[ 

Notice that the presented method can be applied to high order behaviour of per- 
turbation theory not only for the ground state wave function but also for the ground 
state density matrix equaling to the product p = ^(x)^(y) because of the reality of 
the ground state wave function. The density matrix and the k-th order for it can be 
presented in a form analogous to (PS). We can change the variables like eq.(|5]) and 
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find the pk asymptotics as 

k —> oo,x/Vk — > £i,y/Vk — > £ 2 , (13) 

which is given by 

Pk{SiVk,&Vk)~r(k/2)exp 



-k[\(S[^] + S[V>4^]) + l^-\ 



(14) 

where and are the trajectories satisfying the boundary conditions £^(0) = 
Ci)C^(O) — £2- We can also obtain that the trajectories Qi = ^y/X and Q2 = 
are classical, while the parameter A/2 is equal to the sum of the integrals 
fllOD corresponding to these two euclidean solutions. Notice that eq.([L4]) can be also 
obtained by making use of the asymptotics under conditions (|6]), but the derivation 
is not so trivial, because one is in need of estimation of the sum J2n=o ^n(x)^k-n(y) 
under conditions fll3|) . 

Consider the loop-like trajectory Q(r) shown in fig.l. Let Qi(t) = Q(r — 
To), Qz{t) = Q(tq — t),t < 0. This set of trajectories contributes to the asymp- 



totics ( |i~4|) at following values of £i,£2:£i = £2 = Q(0)/y / 25o, while the parameter A 
is equal to 25*o, so that the asymptotics (0) at coinciding arguments takes the form 
T(k/2)So k/2 . Thus, we see that for finding the density matrix large order asymptotics 
at coinciding arguments which is necessary for obtaining the asymptotics of the main 
values of the operators like x ix ...x im ("Green functions") one is in need of making use 
of classical solution Q(r), and the corresponding asymptotics has the form analogous 
to the ground state energy large order behaviour. This result is in agreement with 
the technique developed in H. 

But the results obtained in the presemnt paper allow us to find the asymptotic 
behaviour of the fc-th order of perturbation theory for the main values of the observ- 
ables O being the operators with the kernels 0(x,y,k) depending on k as follows, 
0(x,y,k) = f(x/Vk,y /Vk) exp(kg(x/\/k,y /Vk)). Namely, these quantities being 
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the integrals J dxdypk{x, y)0(x, y, k) can be approximately calculated by making use 
of eq. (|14D and applying the Laplace method. 

One can also calculate other quantities. Namely, the k-th order of perturbation 
theory for the wave functions is presented as a product of a polynomial function 
J2 a ix...i n x% i ■■■ xl n by the Gaussian exponent e~ x2 l 2 . It occurs that the obtained result 
(0) allows us to find the a^...^ asymptotics as h/k, ...,i n /k tend to constants, because 
these coefficients can be written through the contour integrals having the saddle-point 
form under mentioned conditions. 

The investigation of the wave function large order asymptotics gives us a pos- 
sibility to substitute it to perturbative recursive relations, to check the exponential 
approximation (|7]) and to find the pre-exponential factor, as well as the corrections, 
not only by calculating the determinant of fluctuations in the path integral approach, 
but also by direct analysis of recursive relations, which allows us to consider the 
cases of singular quantum corrections |?J to the Hamiltonian and of instanton-anti- 
instanton pair ||. For semiclassical expansion such analysis was made in ||. From 
the recursive relations research one can also confirm the asymptotics ( 12 ) and find 
the function x m it- 

I hope that the suggested method of constructing asymptotics can clarify the 
problem: what classical solutions contribute to the large orders of perturbation theory 
and, in particular, find the correspondence between the existence of this contribution 
and the parity of a number of negative modes around the euclidean solution. The 
correspondence was found in |7j for the ground state energy and confirmed in || for 
the large orders of semiclassical expansion. 

Probably, the suggested method can give rise to a new vision of a multiparticle 
production problem considered recently in many papers, see, for example, [| [II], 
|TT| . For considering this problem it is important to know the behaviour of Feynman 



7 



diagramms with a large number of the external lines and a large number of vertices, i.e. 
the main values of the observables depending on k in the k-th order of perturbation 
theory at larges k which has been already discussed. 

The author is indebted to V.A.Rubakov for helpful discussions. This work is 
supported in part by ISF, grant # MKT000. 
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Figure caption Fig.l. Classical euclidean solution Q(r) and corresponding curve 
in £ -space. 
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